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^ ■ Abstract 

These lectures present a pedagogical introduction to the physics of heavy-flavour pro- 
duction in hadronic collisions. The first lecture gives the theoretical background, with a 
discussion of leading-order calculations and of the effects of next-to-leading-order correc- 
tions. The origin and implications of the large logarithmic corrections appearing at this 
order are presented in an elementary way. The second lecture provides a survey of current 
experimental data on charm and bottom production, and describes their comparison with 
theoretical predictions. We emphasize the role played by some non-perturbative effects in 
the determination of charm distributions, and study the theoretical systematic uncertainties 
which affect our predictions. 
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1 Introduction 

Most of this wonderful School was devoted to illustrating the key role that heavy quarks occupy 
in helping physicists probe the structure of elementary particles and their interactions at the most 
fundamental level. Progress in understanding the main conceptual weaknesses of the Standard 
Model, such as the origin of the family structure or the origin of CP violation, might well come 
in the future from more accurate studies of bottom and top quark decays and interactions. The 
present set of lectures, will address instead less mind-boggling aspects of heavy-quark physics, 
namely the physics involved in their production. Nevertheless, as I hope these lectures will convince 
you of, there is more to learn from a theoretical study of heavy-quark production than just cross- 
sections and predictions for the numbers of b quarks to be expected at the LHC. 

First of all, heavy-quark production provides a benchmark process for the study of perturbative 
QCD. This is because the heavy-quark mass mq defines the scale at which the strong-interaction 
coupling constant a s is evaluated, and mg ^> Aq C d implies that the production properties (at 
least at the inclusive level) should be calculable within perturbation theory (PT). The existence 
of heavy quarks with different masses (m c ~ 1.5 GeV, ~ 5 GeV, mt — 175 GeV), allows 
us to probe perturbative QCD in regions of different Q 2 , where the relative impact of radiative 
corrections and non-perturbative (NP) effects are very different. For example, we know that NP 
effects are very important in turning a c-quark into a charmed hadron (say a D + ). Were we able to 
fully control the charm production within PT, we could then hope to disentangle the NP physics 
from the comparison of PT with data. I shall give some examples of this in the second part of 
these lectures. 

The presence of a heavy quark in the final state is also a flag for very specific production 
mechanisms and selects specific classes of relevant higher-order corrections. This can be useful 
either as a probe of the nucleon structure, or as a test of our understanding of radiative corrections 
in QCD. For example: 

• jy+charm final states probe the strange content of the proton, since the leading production 
mechanism is gs — > Wc. 

• Inclusive b production at high energy probes the gluon density of the proton, since the 
leading process is gg — > bb. 

• Associated production of W and heavy-quark pairs is sensitive to gluon-splitting processes, 
since the dominant production process is the production of a W and an off-shell gluon, which 
then decays to the heavy-quark pair: qq' — > Wg* — > WQQ. 

3 On leave of absence from INFN, Pisa, Italy 



1 



The top quark is the heaviest of the heavy quarks. The history of its discovery and a thorough 
review of its first measured properties have been nicely covered in Bellettini's lectures at this 
School. There you certainly learned how much of a theoretical input will be required in measuring 
important properties such as the top mass or the CKM couplings. The mass, for example, is mea- 
sured in hadronic collisions by reconstructing multi-jet final states from the top decay. Statistical 
errors on mt op achievable at the LHC (<C 1 GeV) are overwhelmed by the systematic uncertainties 
due to the description of top production and decay. In fact a theoretical knowledge of the structure 
of the top-decay jets and their origin is needed to accurately fit the measured mass distributions |1] 
and extract a mass parameter^. A similar need for theoretical input is required when extracting 
the top mass from the threshold behaviour in e + e~ annihilations, one of the goals of future Next 
Linear Colliders |2j . The direct measurement of the CKM matrix element Vtb can be performed by 
detecting production of single quarks in hadron collisions, via the processes qq' — > W* — > tb and 
gq — > tbq 1 (where the exchange of a t-channel W is understood). In both these processes, in fact, 
the production rate is directly proportional to |V^| 2 . Once more, accurate studies of the QCD 
aspect of the production dynamics need to be performed to make a meaningful measurement of 
| Vy, | possible. 

The study of b cross-sections at the Tevatron [[J is helping us to improve the reliability of 
our estimates of b rates at the LHC, where fe's will be used to measure CP violation Q and to 
probe possible 6-meson rare decays, foreseen in several theories beyond the Standard Model, as 
reviewed in this School by Masiero. The expected total cross-sections will define overall rates, the 
kinematical distributions will determine triggering strategies and detector design. 

Heavy quarks are also important signals of possible new physics, in addition to providing large 
backgrounds to searches for new physics: 

Signals: Backgrounds: 

bb 

Wbb 

Wbbjj, tt 
6677 
tt 



H ^bb 


99 


WE -> bb 


qq 


t^b\ + , cx°, tx° 


pp 


b'^jb, Z°b 


pp 


x + xl -> 3 le Pt' s 


pp 



Understanding the characteristics of final states with fe's (and 6-jets) is fundamental in order to 
learn how to isolate signals from backgrounds. 

As a final example of possible uses of heavy-quark production, let me mention the searches 
for quark-gluon plasma done in heavy-ion collisions. Among other probes, final states with J /if) 
and, more generally, DY pairs, will be used to study the properties of matter at high density. 
Charm pairs will provide both signals (cc — > J/if>) and backgrounds (cc — > To make 

robust predictions on the properties of charm produced in high-energy nuclear collisions at future 

4 This is similar to what happens in the case of the Z mass measurement at LEP/SLC, where the Z mass is 
obtained from a fit to the Z lineshape. The lineshape is affected by higher-order QED corrections, mostly initial- 
state radiation. These need to be evaluated theoretically before the Z mass, which is a parameter in the fit, can 
be extracted from a comparison with the data. 
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machines (RHIC and LHC) it becomes important to test our understanding in simpler produc- 
tion environments, such as those provided by current fixed-target hadro- and photo-production 
experiments. Due to the smallness of the charm mass, significant NP corrections are expected. As 
mentioned earlier, this is therefore also an important place to study the interplay of perturbative 
and non-perturbative effects, as will be discussed in the second part of these lectures. 
The structure of the lectures is as follows: 

1. The theory of heavy-quark production: 

• leading order 

• total and differential cross-sections 

• NLO corrections 

2. Phenomenology of heavy-quark production: 

• fixed-target charm production 

• interplay of perturbative and non-perturbative effects 

• bottom production at the Tevatron 

2 Theory of heavy-quark production 

In this Section I will present an elementary introduction to the theory of heavy-quark production 
in hadron-hadron collisions. Production of heavy quarks in e + e~ collisions is to leading-order (LO) 
an electroweak process, and is therefore relatively straightforward to control. The complications 
of QCD arise only at the level of radiative corrections, and become relevant only when one is 
interested in high-precision studies (such as the measurement of Rb at LEP/SLC) or in particular 
aspects of the QCD dynamics (such as the study of fragmentation functions, to which we will return 
in the following). Photoproduction, which is of relevance for several fixed-target experiments and 
for the HERA collider, can be described by a sub-class of the processes appearing in hadronic 
collisions, so we will not consider it in detail here (although we will discuss its phenomenology 
in the second lecture). For a complete and more technical review of the theory of heavy-quark 
production, as well as for a precious list of references to the original papers, see for example 
Nason's work in ref. ||. Additional details and phenomenological applications can be found in 
the excellent book by Ellis, Stirling and Webber ||. 

2.1 Leading order 

Two processes are responsible for heavy-quark hadro-production at the LO in perturbation theory: 

qq -> QQ and gg — > QQ . (1) 
The corresponding diagrams are shown if fig. [I]. Few simple comments are in order: 

• In qq — ■> QQ, the QQ pair is always in a colour-octet state. In gg — ► QQ both colour-singlet 
and octet are allowed. 



3 




a 



sJULSJUU 



JUUUULi 



D 



a 



D 




2 



D3 



Figure 1: Leading- order diagrams for heavy- quark pair production. 

■> The only difference is an overall colour 



The process qq — > QQ is similar to e + e 
factor 

H ^ij^kl E 



1 



(2) 



which when squared and summed over all possible initial and final colour states gives a factor 
C F = (N 2 - 1)/2N (C F = 4/3 for N = 3). 

• The total production cross-section for heavy quarks is finite. In fact no poles can develop 
in the intermediate propagators. This is obvious for the s-channel diagrams, since s > 4m 2 . 
In the case of the t-channel gluon exchange (say in diagram D\ of the figure), labeling the 
lower-line gluon and quark momenta as p\ and Q, respectively, one obtains in the CM frame: 



(pi - Qf -m 2 = -2 Pl ■ Q 



- S -(l 



[3 cos 9) 



(3) 



with P the heavy-quark velocity (f3 = \Jl — ^f- ) and 9 the scattering angle. Therefore the 
propagator never vanishes: 



2 Pl ■ Q > £(1 -/3) = — J— -(1 - /? 2 ) = > m 2 



(4) 



Therefore m 2 is the minimum virtuality exchanged in the t-channel, and the total cross-section 
is finite at LO in PT. This would not be the case for massless quarks. The minimum transverse 
momentum transferred in the t-channel sets the scale for the argument of the strong coupling 
constant a s . Since this scale is of the order of the heavy-quark mass, and since this is (by 
definition of heavy quark) much large than Aqcd, perturbative calculations should be reliable. 
Once again, this could not be the case for light quarks: the total production cross-section for up 
or down quarks in not calculable in PT! 
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The evaluation of the matrix elements can be done using standard Feynman diagram tech- 
niques. The results are: 



1e 

1 „ 



M{qq - QQ) 
M(gg - QQ) 



V 2 2 Pn 

2^^! +r 2 +- 



/ 1/ 



2VJV Vrir 2 



-2iV' 



Ti + Tn 



P 



4riT 2 , 



(5) 
(6) 



where: V = N 2 — 1 is the dimension of the SU(N) gauge group, i.e. the number of gluons, and 

_ Pi 2 ■ Q 1 =F P cos 4m 2 



Tl, 2 



(7) 



Exercise 1: show that in the limit TiT 2 — > (which can be achieved at very high energy when 
(3 — > 1) the gg — > QQ matrix element squared goes at most like 1/(tit 2 ), and not like 1/(tit 2 ) 2 , 
as the eq. @ naively indicates. 

Exercise 2: show that the term inside the first parenthesis of eq. (pi) is always positive. 



The total partonic cross-sections are obtained by integrating over the two-body phase space: 

"*» s h wm> ^ - p -» = 1 (s) 2 " < 8 > 

giving the following results: 

- + (9) 

Hgg^QQ) = § (^) ^p [*m (p 2 + 2^( P + 1 

+2(V - 2)(1 + p) + p(6p - iV 2 )l l£(/3) (10) 



where 



Comments: 



At large s the gg rate vanishes more quickly. 

Notice that the threshold suppression in the gg — > QQ process goes like 

m + P) = ^~P (2 + P) * 2(1 - t )(1 + P - ) = 2(1 - 4^) (12) 

As a result threshold effects vanish very quickly as soon as s > Am 2 . This suppressed 
suppression is related to the spin- 1/2 of quarks. For scalar particles (e.g., squarks) the 
phase-space suppression is stronger (~ /3 3 ). 
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This is important for top production, which for the measured value of the top mass is dominated 
at the Tevatron by qq annihilation close to the kinematic threshold. Including the reduction in 
rate due to the fewer degrees of freedom of a scalar particle relative to a spin-1/2 one, the squark 
cross-section turns out to be approximately 1/10 of that of a quark with the same mass, if the 
dominant production mechanism is qq annihilation. 



Examples: 

• The relative productions rates for heavy quarks of mass mi and m 2 behave, at high energy, 
like: _ 

v(99 -> QiQi) ™ l _ bgjm|/m|) 
H99 Q2Q2) log(s/m!) 

For example, 

a{ " - ^ ~ 0.7 at v^=100GeV 



a(gg -> cc) 
For the qq process we have instead: 

£(gg -» Q1Q1) 

a(qq -> Q 2 Q 2 ) 

with 



(14) 



1 - 0{mi/s 2 ) (15) 



- bb \ ~ 0.99 at v / I=100GeV 
<r{qq — > cc) 



(16) 



In general, we are interested in hadronic rather than partonic rates. The obtain these we need to 
convolute the partonic cross-sections with the parton densities in the hadron. It is then useful to 
parametrize the final state in terms of the transverse momentum (p±) and of the rapidities (y) of 
the produced quarks. The resulting phase-space, including the integration over the momentum 
fractions of the hadrons carried by the initial-state partons x±^, is then: 

dfa = ^ dx 1 dx 2 -^^ s ^|=o (2*-) V(P in - P out ) (17) 



Use: 



dx x dx 2 5(E- m - E out )5(P-: Q - P: ut ) = — ^— = (18) 

Z - C/ bcam D had 



dQ z 1 Q° + Q z 9 

— = dy, y = - log ^ Q - ^ z (-> 77 = - log(tan -) for m = 0) (19) 



to get: 



77x1X2 / 1 



2 



dydydpr (20) 



It is easy to rewrite s in terms of y, y and of the transverse mass m\ = p\ + m 2 using: 

<5o = Tfi r p cosh?/ 
Q 2 = m T sinhy 



(21) 
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Figure 2: Total production cross-sections for charm, bottom and top quark pairs, in pp 
collisions. 

This gives: s = 2m j, {1 + cosh(y — y)} and finally: 

W T = £f [i + coshfa-,)]^ x E««>WiWIW« - «)l 2 < 22 > 

For a fixed value of pt (and therefore of m^), the rate is heavily suppressed when \y — y\ becomes 
large. Therefore Q and Q tend to be produced with the same rapidity. This has important 
consequences for the design of detectors aimed at collecting large numbers of bb pairs, as required 
to perform CP studies. In fact, it implies that once a detector has registered the b quark, there is 
a large probability that the b will be nearby in rapidity. This is fundamental for detectors covering 
the forward region only, since it shows that there is no need to build a backward spectrometer as 
well to efficiently collect the heavy-quark pair. The gain in acceptance for the pair would only be 
a factor of 2, and not 4 as naively expected in absence of such strong rapidity correlations! 

2.2 Some results 

We discuss here some simple applications of the cross-section formulae obtained in the previous 
section. They can be convoluted numerically with standard parametrizations of the parton densi- 
ties inside the proton, to obtain hadronic cross-setions at different energies and for different quark 
masses. In fig. ^| I show the total production cross-sections in pp collisions, as a function of the 
hadronic center of mass energy y~S, for charm, bottom and top quark pairs. Notice the different 
units used in the case of the top quark relative to those used for charm and bottom! Top cross 
sections are suppressed by several orders of magnitude even at the highest energies reachable in 
the foreseeable future. When the partonic CM energy grows, however, the suppression due to 
the large mass becomes less and less important, as discussed in the previous section. This can 
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pt (GeV) 




Figure 4: Inclusive pr distributions for charm, bottom and top quarks pairs, in pp 
collisions at y/S = 20 TeV. 

be seen in fig. [3], which shows the transverse momentum distributions at v^S* =1.8 TeV. Notice 
that the charm and bottom rates become almost equal as soon as pt > 40 GeV (which means 
y/§ > 80 GeV). In the case of the top, we need to go out to p? > 400 GeV. It is instructive to see 
how this picture changes at higher hadronic energies. This is shown in fig. | for VS = 20 TeV. 
Notice that at pt = 400 GeV there is still a large difference between the c, b rates and the top rate. 
This is because, contrary to the case of 1.8 TeV, at the LHC the production rate is dominated by 
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gg annihilation. As we showed earlier, the mass effects vanish more slowly for gg processes than 
for the qq annihilation process. More details on the comparison of the predicted production rates 
with the available experimental data will be given in the second lecture. For the time being, we 
shall move on to the discussion of the corrections beyond the Born level. 



2.3 Next-to-leading order corrections 

Next-to- leading-order (NLO) corrections come from two sources of 0(o? s ) diagrams: 



Real Emission Diagrams 



OJUUUUU 



JUUUUUU 




Virtual Emission Diagrams 






In the first case, the corrections come from the square of the real emission matrix elements. In 
the second case, from the interference of the virtual matrix elements (of 0(g A )) with the Born 
level ones (of 0(g 2 )). Ultraviolet divergences in the virtual diagrams are removed by the renor- 
malization process. Infrared and collinear divergences which appear both in the virtual diagrams 
and in the integration over the emitted parton in the real-emission processes cancel each other, or 
are absorbed in the initial-state parton densities. The complete calculations of NLO corrections 
to the production of heavy-quark pairs in hadro- and in photo-production were done in ref. [7J § 
(total hadroproduction cross-sections), ref. |], [ITJ (one-particle inclusive distributions in hadropro- 
duction), ref. [[□], [L2| (total and one-particle inclusive distributions in photoproduction), ref. |p~3| 
(two-particle inclusive distributions in hadroproduction) and ref. [13 (two-particle inclusive dis- 
tributions in photoproduction). The explicit results are too complicated for me to reproduce them 
here. So I will limit myself to a qualitative discussion only. 

NLO corrections should improve our knowledge of the production cross-section. In particular, 
uncertainties related to the renormalization (/%) and factorization (fip) scale dependence should 
be reduced^. There is evidence however that the NLO is not sufficient to get accurate estimates, 
since a large scale dependence is still present. This is shown is fig. |5], which shows the scale 
dependence of the inclusive pt distribution of b quarks at the Tevatron. Large scale dependence 
is a symptom of large NNLO corrections A full NNLO calculation goes beyond our ability. So, we 
are stuck. 



5 The renormalization scale /zr is the energy scale used in the evaluation of a s 
scale used in the evolution of the parton densities. 



The factorization scale /zf is the 
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Figure 5: Scale dependence of the inclusive pt distributions for bottom quarks, in pp 
collisions at \f~S~ = 1.8 TeV. /xr = fip = fi. 

We can however make some progress in this direction by trying to understand what are the 
sources of these large NNLO corrections. We can classify them as follows: 

• "Accidentals" (e.g. large numerical coefficients, 7r 2 's, etc., which should not affect the scale 
dependence) . 

• Large logarithms due to the coexistence of different scales. E.g.: 



TXT 



, Pt 
log — 



The origin of these logarithms is known, and will be discussed in the remaining of this lecture. 
They are universal, namely they are associated to a universal behaviour of the cross-section in 
some particular kinematic configurations. Because of this, their structure can be predicted at all 
orders of PT, and in principle they can be resummed. One then hopes that these resummed large 
logarithms capture all the essential features of the higher-order corrections, and help to improve 
the behaviour of the perturbative predictions, without the need of a complete NNLO (or even 
higher-order) calculation. 

Before we start analysing the sources of these logarithms in detail, some general comments are 
in order: 

i) No collinear singularities appear when gluons are emitted from the final state heavy quarks, 
since they are screened by the quark mass. Therefore, contrary to the case of a light parton, 
da/dpr for a heavy quark is a well-defined quantity in NLO. For light partons, one would 
encounter a collinear singularity and would have to introduce a fragmentation function, 
usually not calculable from first principles. 



10 



At large pr, nevertheless, large log(pr/m) factors appear, signalling the increased probability 
of collinear gluon emission. At large p?, the massive quark looks in fact more and more like 
a massless particle. These logarithms can be resummed using the fragmentation function 
formalism. Contrary to the case of light quarks, however, the fragmentation function for a 
heavy quark can be calculated from first principles [113], |TB |. 



iii) New processes appear at NLO which drastically change either the s dependence of the 
cross-sections, or the kinematical distributions. 

In the following few subsections I will present some examples of where these large corrections 
arise, and how they can be accounted for. The discussion will be rather qualitative, and is just 
meant to give an idea of the main issues and to introduce a nomenclature which is often used 
in the literature and heard in seminars. I hope these examples will help you demistifying some 
possibly esoteric language. For more detailed treatments, look at the literature quoted. 



2.3.1 i-channel gluon exchange diagrams 

Diagrams with a gluon exchanged in the t-channel alter the high-energy behaviour of the total 
cross-section [0. Consider for example the following diagram: 




Let us call s the incoming total energy squared, and s' the energy squared of the hard subprocess 
involving gluon 2 and 1'. Let us concentrate on the phase-space region where s' <C s. Since the 
probability to find the gluon 1' with momentum fraction x inside gluon 1 is proportional to a s /x, 
the approximate contribution from this process to the heavy-quark cross-section is 



Using 



VI- 



Therefore, 



^oox ~ / „ , <&- *(1'2 - QQ) (23) 



4m 2 /s X 



a 2 a 2 
s 

s' xs 



1 dx a 2 al 



^NLo(gg^QQg) «->«> s 

a s~i ( 24 ) 



VLo(gg -> QQ) ™ 2 



n 



At very high partonic energies, the NLO process will dominate over the LO one! Of course 
production at large energy is suppressed in the total hadronic cross-section by the convolution 
with the steeply falling parton densities. To quantify the impact of these effects, we can take for 
example a gluon density inside the proton with the form: 



/(*) 



A 

x l+8 



with 5 < 1 It is then easy to show (Exercise!) that 

o"nlo(pp -> QQ) 



«slo gi ^ if51o g3 ^«l 



*to(RP-QQ) [ as i±s if5 i og _^»i 



Therefore, for small 5 the NLO correction becomes larger than the LO term as soon as a s log S/ 4m 2 > 1. 
At the Tevatron (and for m = rrib ~ 4.5 GeV), this number is indeed of order 1! 

These large logarithms are known as "log -" terms, where x ~ JAm 2 /S is the average value of 



the hadron energy fraction needed to produce the pair. These logl/x terms appear at all orders 
of PT, but fortunately techniques are known to resum them W7\, [TBI]. Notice, however, that they 



would not be present if the gluon density were steeper than 1/x. For example, 5 = 0.5 would 
give only a moderate NLO contribution. The apparent paradox of this different behaviour for 
different values of 5 is resolved by the observation that these large log 1/x terms can be resummed 
in a universal fashion and absorbed in an improved evolution equation for the gluon density (the 
BFKL evolution equation). As a result of this evolution, a gluon density with 5 = immediately 
develops a 5 > 0. One can therefore work with a BFKL resummed gluon density, and forget about 
large log(l/x) corrections to the heavy-quark production cross-section. 

More accurate estimates of the effect of log(l/x) resummation for b productionat the Tevatron 
can be found in |p.8|] , where it was concluded that corrections of up to 30% in addition to the NLO 
calculation can be expected. 



2.3.2 Gluon-splitting contributions 

The so-called gluon-splitting contributions come from diagrams of the following type: 




In this case we should focus on the region of phase-space where the t-channel gluon is very hard and 
far off-shell, while the splitting gluon has a virtuality much smaller than its transverse momentum. 
These diagrams then give rise to potentially large logarithms. We can use a fragmentation-function 
language to estimate the contribution of these processes to heavy-quark production. One factorizes 
the process into a gg — > gg (or qg — > qg) scattering, with one final-state gluon off-shell, followed by 
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the evolution of the off-shell gluon into a QQ pair. Higher-order corrections, such as the emission 
of multiple gluons from the final state gluon, can be included by convoluting the light-parton 
cross-section with the g — > Q fragmentation function. This approach is similar to what done in 
the case of inclusive light-hadron production inside a jet, with the difference that the g — > Q 
fragmentation function can be calculated from first principles within PT |]16|| . This will account 
both for the increase in probability that the gluon will turn into a heavy-quark pair (due to the 
larger number of gluons present when one considers the multi-gluon emission), and for the slow- 
down of the heavy-quark momentum due to the loss of energy in the multi-gluon radiation process. 
The fragmentation function for a heavy quark inside a gluon jet with transverse energy p± obeys 
the evolution equation: 

dD(x,q 2 ) a s f 1 dz 2 2 x y 

where G(z, q 2 ,pT 2 ) is the density of gluons of virtuality q 2 carrying a momentum fraction z inside 
a gluon of maximal virtuality pt 2 , and Pq 9 (x) is the Altarelli-Parisi kernel given by 

The heavy-quark multiplicity (which, for small average multiplicity, coincides with the probability 
to find a quark in the gluon jet) is given by the first moment of the fragmentation function^: 

N= [ dzD(z) (26) 

which obeys the equation: 



-G(z,q 2 ,p T 2 )P Qg (-) (25) 



o 



dN _a s (1) 



d log q 2 2tt 

with G(i) and Pqg first moments of G(x) and Pq 9 (x). G^ ) {q 2 ,pr 2 ) is the number of gluons with 



virtuality q 2 in the jet, and 



z) 2 ) 



1 

3 



At the first order, neglecting the g 2 -dependence of the gluon multiplicity and therefore assuming 
G(i) = 1, we get 

HQQ) ~^log(4) (27) 
D7T \m 2 J 

where p\ and m 2 are the maximum and minimum gluon virtuality. In the case of light quarks, 
m = and this calculation would not be possible. A collinear divergence would develop, associated 
to the vanishing minimum gluon virtuality, and one would need to introduce a non-perturbative 
fragmentation function. Once more, this shows the advantage of working with heavy quarks! 

The logarithm in eq. (|27]) signals the growth in probability to find a heavy quark inside a 
large-pr gluon. It is yet one more potentially large logarithm that was not present at leading 
order, and which appears for the first time at NLO. Similar large logarithms will appear at higher 
orders of perturbation theory (associated to processes where more gluons are radiated in the final 
state). These logarithms are resummed by solving exactly the evolution equation for the g —>■ Q 
fragmentation function, eq. 



6 The n-th moment of a function f(x) is defined by dx x n 1 f{x). 
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2.3.3 Flavour-excitation processes 



An analogous class of NLO corrections is given by the so-called flavour- excitation diagrams. These 
can be thought of as initial-state gluon-splitting processes: 




The relevant region of phase-space in this case is the one with the heavy-quark propagator close 
to the mass shell. These processes are analogous to the singlet contribution to F2 in DIS. It is 
known that these effects can be reabsorbed in the AP evolution of the parton densities fli~9 |. In 
our case, this is equivalent to defining a heavy-quark density inside a proton by: 



dD Q (x,q 2 ) 
d log q 2 

D Q (x,m 2 



§/ifG(,,g 2 )P Q9 (f) 



(2* 



0. 



It is easy to get a crude estimate of the "heavy-quark density" by assuming a simple functional 
form for G(x), for example, G(x) = A/x. Then: 



d log q 2 



27T Jx Z \Z 



p, 



27T Jx 

a. 



1 dz A 1 



z x/z 2 



QgK 



z 2 + (l 



6tt 



G{x) 



Assuming G(x, q 2 ) slowly varying in the q 2 > m 2 region: 




Since the NLO corrections contain these "structure function-like" piece, it would be incorrect to 
add to LO processes contributions from gQ — > gQ. In fact, this would imply a double-counting 
of the NLO flavour-excitation diagrams. Of course, the use of the gQ — > gQ process allows 
one to capture higher-order effects included in the AP evolution equations, which are usually 
solved to generate the heavy-quark densities found in standard Fortran libraries of parton density 
functions. The use of the NLO flavour-excitation diagrams reproduces instead more faithfully 
the exact kinematics and correlations of the flavour-creation process in the region next to the 
threshold. 

The example I will give here describes the comparison between the two approaches when 
applied to the simpler case of associated production of a 7 and heavy quarks. In this case I 



14 



compare the results of the two calculations of the photon p± spectrum obtained by using either 
of the these two processes: 



Born 




k/V\/V 



The calculation done using the LO approach, which uses as an input the b structure function 
taken from a given set of parton distribution functions (CTEQ1M, in this case), is shown in 
fig. H as a solid line. The calculation done by evaluating the 0{a s 2 a-EM) matrix element is shown 
by the dashed line. The ratio of the two distributions is given in fig. 0. Notice that the two 
calculations agree remarkably well down to low values of px, indicating that even the structure 
function approach can be used to reproduce rather accurately the onset of the mass threshold 
(this is because we are looking at a rather inclusive quantity. Were we looking at more exclusive 
quantities, or at correlations between the b and the b in the final state, the LO approach would 
not even be usable: in this approach the antiquark is always integrated over and we cannot keep 
track of where it went!) At large px some differences between the two calculations emerge, due to 
the absence of higher-order terms in the NLO calculation. These higher-order terms are included 
in the structure-function approach, since the b structure function itself has been evolved to large 
Pt using the exact Altarelli-Parisi evolution equation. The ideal calculation, therefore, would 
merge the two approaches by smoothly interpolating betwen the two regimes. This was done in 



refs. [20, ^TJ, and the results will be discussed in the next lecture. 



2.4 Heavy-quark fragmentation 

At the end of the perturbative evolution, one still expects hadronization effects to alter the spec- 
trum of the heavy quark while it transforms into a hadron. These effects are usually described 
in terms of a non-perturbative fragmentation function. The fundations of this approach lie in a 
factorization theorem: 

dP T J z dp t 

where Pt is the hadron momentum, and da/dp t includes the perturbative part of the heavy-quark 
fragmentation function, namely the resummation of the log (pr/m) terms The factorization 
theorem ensures that, up to corrections of order A./ Pt, with A a scale typical of the hadronization 
phenomena, the function F(z) for a given heavy flavour is universal. Although non-perturbative 
in nature, some general features of fragmentation functions can be extracted from first principles. 
For example, one can establish a relation between the effects of hadronization on heavy quarks of 
different masses: 

l-x B =(—)(l-x D ) (29) 
\m b J 
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Figure 6: The px distributions for central photons produced in asociation with a b quark, 
in pp collisions at \fS = 1.8 TeV. The results of the two different calculations discussed 
in the text are displayed. 
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Figure 7: The ratio between the two curves shown in the previous figure. 
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where x is any parameter such as (x) (the average value of x) or x max (the value of x at which 
the fragmentation function peaks). This follows easily from phase-space considerations. One can 
idealize the fragmentation process with a transition of a slightly off-shell heavy quark (possibly in 
the field of nearby soft gluons) decaying to a heavy hadron H plus a light system II (e.g. pions): 



Q 




In this case: 

P H : P n = z H : z u = M H : M n (30) 

and, clS du result, 

1 — z ~ 1 ~ 31 

M H + M n M H y ! 

A typical example of parametrizations for non-perturbative fragmentation functions is the Peter- 



son, Schlatter, Schmitt and Zerwas form f22| : 



dN 1 

— — oc 7j (32) 

dz x [l - l/x-e/(l -x)] 2 

Here e is the only parameter describing the non-perturbative part, and scales as follows: 

9 

TT) 

e ~ -f- (1 - 2) ~ Vi (33) 

7T7,n 



Another interesting form was proposed by Nason and Colangelo |23 

^ oc (1 - ^) Q / (34) 

In this case the relation between the values of a and /3 for different values of the heavy-quark 
mass (say charm and bottom hadrons) can be extracted by imposing the two relations: 

l-(x B ) = ^(l-(x D » (35) 
m b 

l — x B = —(l-xjj), (36) 

which can be easily solved analitically (Exercise!) in terms of a and j3. The values of the pa- 
rameters characterizing the non-perturbative fragmentation functions can be extracted by fitting 
some set of data. The factorization theorem will then allow us to use the same parameters to 
predict cross-sections for other data. The best place to extract information on heavy-quark frag- 
mentation is e + e~ collisions, since there at least the energy of the primary heavy quark before the 
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Figure 8: Fits to the charm fragmentation function measured at ARGUS (left) and 
OPAL (right). 

perturbative and non-perturbative fragmentation processes occur is well determined, being equal 
to the beam energy. Recent fits to charm production data from the old ARGUS experiment p4j 
at DESY (y/s = 10.6 GeV) and the more recent OPAL experiment [f25| at LEP (y/s = 91.2 GeV) 



have been performed in ref. ||26|| . The results are shown in fig. || The parameters extracted from 
the fits to the data of the two experiments, using different values of Aqcd, are contained in the 
following tables (for ARGUS and OPAL, respectively): 



A 5 =200 MeV 300 MeV 
e c =0.019 0.011 
a = 1 
= 6.7 



A 5 =200 MeV 300 MeV 

e c =0.015 0.008 
a = 0.9 
P = 6.4 



Few comments should be made: 

• There is a very good agreement between the value of the parameters extracted from the two 
experiments. This confirms the validity of the factorization theorem. 

• The value of the parameters depends strongly on the input value of Aq C d- This is reasonable, 
since the perturbative part of the heavy-quark evolution will depend on the strength of the 
perturbative coupling constant a s . A larger value of a s implies a larger amount of energy 
radiated off in the form of gluons. This leaves less energy to the quark before it hadronizes, 
and calls for a harder non-perturbative fragmentation to reproduce the data than if a small 
value of a s were used to start with. This correlation is faithfully reproduced by the fits shown 
above. Notice also that the value of the Peterson parameter e c which is usually obtained 
by using the LO result for the perturbative part of the fragmentation function (e c = 0.06, 
|2"6"| , |2"Tf), is much larger than what extracted from these NLO fits (e c = 0.01 — 0.02). 

The relations given in eqs. (|3~5l - j36l ) allow us to predict the value of the parameters governing 
the non-perturbative fragmentation of b quarks. When this is done, one can predict the 6-hadron 
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Figure 9: Comparison of OPAL b- fragmentation function data with fits to different 
fragmentation function forms (upper figure) , and with predictions based on the convolution 
of the b-quark NLO fragmentation function and the Colangelo-Nason non-perturbative 
fragmentation function, with parameters extrapolated from fits to the charm fragmentation 
function (lower figure). 



fragmentation function at LEP. The comparison of data |28| and theory (using the Colangelo and 
Nason fragmentation model, and the a and (3 values extracted from the charm fits) is shown in 
the bottom plot of fig. |9|. Excellent agreement is seen, again a good test that we are on the right 
track! 



3 Phenomenology of heavy-quark production 

The ultimate test of our ability to describe the heavy-quark production processes rests on the 
comparison with actual data. In the remaining part of these lectures I will present a survey of 
current experimental results, and how they compare to theory. Unless otherwise indicated, theory 
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will refer to NLO calculations. A more complete review can be found in ref. 
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3.1 Fixed target production 

Heavy-flavour production has been studied extensively in fixed-target experiments, with both 
hadron and photon beams. The typical centre-of-mass energy is in the range 10-40 GeV, where 
the bottom cross-section is rather small. Therefore most of the available data are on charmed- 
hadron production. Total cross-sections, single-inclusive distributions, correlations between the 
quark and the antiquark have been measured in both hadro- and photoproduction. The theoretical 
apparatus of perturbative QCD is in this case at its very limit of applicability, because the charm 
mass is very close to typical hadronic scales. Thus, effects of non-perturbative origin will very 
likely play an important role. Conversely, it is hoped that these effects may be better understood 
by studying charm production. Although modern fixed-target experiments have considerably 
improved the situation, many open problems remain in this field. All experimental results are 
in qualitative agreement with perturbative QCD calculations, thus supporting the "hard" nature 
of charm-production phenomena. However, several quantitative deviations from pure QCD are 
observed. It is interesting to see whether simple models of non-perturbative phenomena, such 
as fragmentation effects and intrinsic transverse momenta, may be sufficient to explain these 
deviations. I will discuss these problems at length. I will not discuss, instead, the subject of 



^-asymmetries. A discussion of this issue can be found in ref. [29], and recent relevant data (and 
references to previous literature), in ref. p0 . 




Figure 10: Pair cross-sections for b and c production in ir N collisions versus experi- 
mental results. 
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3.1.1 Total cross-sections 



In fig. pXj I plot the cc and bb cross-sections, computed in QCD at NLO, as functions of the beam 



energy, for n~N collisions. The same quantities are shown in fig. |TT] for a proton beam. The 
cross-sections are calculated using the parton distribution sets of ref. [31] for the nucleon, and 
the central set SMRS2 |E| for the pion. The default values of the charm and bottom masses 
are 1.5 and 4.75 GeV respectively, and the default choices for the factorization scale /ip and the 
renormalization scale are 

fi F = 2m c , fi R = m c and fi F = fi R = m b (37) 

for charm and for bottom. 

The bands in the figures are obtained as follows. We varied [ir between half and twice the 
central value. The factorization scale /xf was also varied between rrib/2 and 2m{, in the case 
of bottom, while it was kept fixed at 2m c in the case of charm. This is because the adopted 
parametrizations of parton densities are given for Q 2 larger than 5 GeV 2 . The bands shown in 
the figures are therefore only an underestimate of the uncertainties involved in the computation 
of charm production cross-sections. Considering independent variations for the factorization and 
renormalization scales does not lead to a wider range in the bottom cross-section for the energies 
shown in the figures. We also show the effect of varying m c between 1.2 GeV and 1.8 GeV, and 
rrib between 4.5 and 5 GeV. 




Figure 11: Pair cross-sections for b and c production in pN collisions versus experi- 
mental results. 



The proton parton densities of ref. |3T]] are available for a wide range of A QCD values, cor- 
responding to a s (m z ) values between 0.105 and 0.130. The bands shown in figs. [TO and [TTJ for 
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bottom production are obtained by letting A QCD vary in this range. In the case of charm, values 
of A QCD corresponding to a s {m z ) above 0.115 induce values of a s {m c ) too large to be used in a 
perturbative expansion. For this reason, the upper bounds on charm production cross-sections 
are obtained with a s {m z ) = 0.115. We point out that, by varying A QCD , one is forced to neglect 
the correlation between A QCD and the pion parton densities, which were fitted in ref. [S2| with 
Af = 122 MeV. 

As can be seen, experimental results on total charm cross-sections [^, ^B], [J7], |3S| are 

in reasonable agreement with theoretical expectations, if the large theoretical uncertainties are 
taken into proper account. We can see that the hadroproduction data are compatible with a value 
of 1.5 GeV for the charm-quark mass. In the case of bottom production f39|| , the spread of the 
experimental data is almost as large as that of the theoretical predictions. The two results in 
proton-production [40, 41] differ at the level of 2.5<r. This is unfortunate, since the beam energy 
(E p = 800 GeV) is almost the same of the forthcoming HERAB experiments. A more consistent 
experimental estimate of the production cross-section would have provided more reliable estimates 
of the b rate at HERAB. 




Figure 12: Pair cross-sections for c production in jN collisions versus experimental 
results. 



Total cross-sections for charm production have also been measured in photoproduction experi- 
ments. In fig. [T2| the relevant experimental results of refs. [|2[ |33], [0], are shown in comparison 
with NLO QCD predictions. As can be seen, the theoretical uncertainties are smaller in this case 
than in the hadroproduction case. Again, a charm mass of 1.5 GeV is compatible with photopro- 
duction data. It should be stressed however that some of the experimental results are incompatible 
with one another. Until these discrepancies are resolved, it will not be possible to use the data to 
constrain physical parameters. For example, while the E687 data are inconsistent with a charm 
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Figure 13: Experimental p% distribution compared to the NLO QCD predictions (m c = 
1.5 GeV), with and without the inclusion of non-perturbative effects, in jN collisions. 

mass of 1.8 GeV, this mass value cannot be excluded because of the E691 data. 
3.1.2 Single-inclusive distributions 

While total cross-section measurements are certainly important, it is equally important to explore 
in more detail the production properties of heavy quarks. In addition to the theoretical interest, 
this is relevant to guarantee that some of the experimental systematics (such as the detector 
acceptances and efficiencies) are under control. When experiments measure total cross-sections, 
they often do so by probing a fraction of the full available phase-space, and correcting later for 
the remaining part of phase-space. This is done by using a theoretical modeling of the differential 
distributions. Checking that these agree with the data is therefore a fundamental test. 

Many experiments have measured single-inclusive x F and pt distributions for charmed hadrons 



in nN collisions 0, [46], ^7], |48|, and in pN collisions |36], |37|, |38|, |48|| . Distributions are expected 
to be more affected by non-perturbative phenomena than total cross-sections. For example, an 
intrinsic transverse momentum of the incoming partons, and the hadronization of the produced 
charm quarks, may play an important role in this case. We will therefore try to assess the impact 
of such phenomena by means of simple models. 

I start by presenting in fig. [13] the p\ distributions measured in photon-nucleon collision by the 



E687 E3| and by the E691 |^3], 44 1 collaborations. I also show the NLO QCD prediction for bare 
quarks (the solid histogram), which as you can see is significantly harder than the dataj]. This is 
reasonable, since the non-perturbative fragmentation of the charm quark into a charmed hadron 
will reduce the hadron momentum relative to the quark one. As discussed in the previous lecture, 
we can describe this slow-down by including a non-perturbative fragmentation function. The 
dashed and dotted histograms in the figures include the convolution with a Peterson fragmentation 
function (as well as with an intrinsic transverse momentum, k T , as will be discussed in a moment). 



We chose here the standard value of e c = 0.06, which was used in p9[ . The agreement with the 
data is good. 



7 In this subsection the relative normalization of experimental distributions and theoretical curves has been fixed 
in order to give the same total rate. 
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Figure 14: The single-inclusive p\ distribution measured by WA92, compared to the 
NLO QCD predictions, with and without the inclusion of non-perturbative effects. 

Let us consider now the case of hadropro duct ion, and focus on a recent measurement of the 



single-inclusive differential cross-sections performed at CERN by the WA92 collaboration [[35 
which uses a ir~ beam of 350 GeV colliding with isosinglet nuclei. Consistent results are obtained 
by using the E769 collaboration pi data, using pion, proton and kaon beams of 250 GeV on 



isosinglet targets, and more recent E706 data |5T] using 515 GeV pion beams. 

In fig. |TJ] I show the comparison between the single-inclusive p 2 distributions measured by 
WA92 and the theoretical predictions. The solid curve represents the pure NLO QCD prediction 
for charm quarks, and as you see it almost exactly goes through the experimental points. But this 
is not necessarily a good sign, since we just learned in the study of photoproduction that we need to 
include an important effect due to the non-perturbative fragmentation. The factorization theorem 
discussed in the previous lecture tells us that these effects should be present regardless of the 
nature of the beam, and once we find out that they must be there in photoproduction, we should 
include them in hadroproduction as well. We would therefore expect that the final theoretical 
prediction, including the fragmentation effects, will not agree with the data. Fortunately there is 
one more effect of non-perturbative origin which comes to our rescue. That is the famous intrinsic 
k T . Partons inside the proton cannot have rigorously zero transverse momentum, because of the 
Fermi motion. It has been known for a long time, mainly from the measurement of Drell-Yan 
events, that quarks in the proton have an average intrinsic k T of the order of 700 MeV. The dot- 
dashed line in fig. |TJ] is obtained by convoluting the quark-level results with the same Peterson 
fragmentation function used in the case of photoproduction, and with a gaussianly distributed 
intrinsic k T , with (k^) = 1 GeV 2 . This is a reasonable choice for (fc 2 ), since gluons can be 
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expected to have a braoder k T distribution than quarks. 

As you can see the effects of fragmentation and of the intrinsic k T almost completely cancel 
each other, and leave a p\ distribution which is in fair agreement with the data. If we use the 
value of the Peterson parameter e c extracted from the recent fits described in the first lecture, 
e c = 0.01, we obtain an even better agreement (see the dotted line in fig. [TJ|). Notice that in the 
photoproduction case the p\ distribution is less sensitive to the choice of (A^). This is because in 
photoproduction only one of the two initial-state particles (i.e. the gluon) can have an intrinsic k T , 



while in the hadropro duct ion case both of them do. This reduced sensitivity is shown in fig. [L3 
the curves with (k%.) — 2 GeV 2 and (A; 2 ) = 1 GeV 2 are not significantly different from one another. 

Having determined the values (or at least a range of values) for the non-perturbative parameters 
which we expect should describe the production of charm quarks at large Pt, we must now proceed 
and see whether the remaining possible distributions that can be predicted do agree with the data 
or not. If they did not, this would be an indication that there are extra phenomena at play which 
we have not accounted for. 

3.1.3 Correlations in fixed-target production 

A distribution which is very sensitive to the presence on a non-perturbative fc T -kick is the azimuthal 
correlation between the quark and the antiquark. This is clear because already at LO, where the 
quarks are produced back-to-back in the transverse plane and the A0 distribution is a sharp peak at 
180°, any transverse kick to the system will drastically change the topology of the final state. Many 
experimental results on correlations between charmed particles in hadro- and photoproduction 
have been obtained by different experiments (see for example refs. |p2| , |53| , |38| , [54] , |55| , [49], |5f 



these reported distributions of the azimuthal distance between the charmed hadrons, the rapidity 
difference, the invariant mass and the transverse momentum of the pair. A detailed comparison of 
these results with QCD predictions is performed in ref. P5|. More recently, new measurements of 
the azimuthal distance and pair transverse momentum for charmed mesons have been presented 
by the WA92 collaboration ||58|1 . In what follows we will focus on the distribution of A<j), defined 
as the angle between the projections of the momenta of the pair onto the transverse plane, and 
of the transverse momentum of the pair p T (QQ). We will discuss whether NLO QCD predictions 
can describe the available experimental data. 

In leading-order QCD the heavy-quark pair is produced in the back-to-back configuration, 
corresponding to A0 = tt and p T (QQ) = 0. NLO corrections, as well as non-perturbative effects, 
can cause a broadening of these distributions, as illustrated in ref. |59|] and |29| . 



We have chosen, as an illustration for hadroproduction, the cases of the WA75 and the WA92 
results, which have both been obtained in it~N collisions at the same energy, = 350 GeV. Let us 
first consider the A0 distribution. In fig. [31| we show (solid curves) the NLO result superimposed 
on the data of the two experiments. The charm mass was set to its default value, m c — 1.5 GeV. 
In both cases, we see that the experimental data favour a much broader distribution than the pure 
NLO QCD result for charm quarks. 

One should, however, take into account also non-perturbative effects, as in the case of single- 
inclusive distributions. The dashed and dotted curves in fig. [T5] correspond to the NLO prediction, 
supplemented with the effect of an intrinsic transverse momentum with (A; 2 ) = 0.5 GeV 2 and 
(/c 2 ) = 1 GeV 2 , respectively. We see that with (/c 2 ) = 0.5 GeV 2 it is impossible to describe the 
WA75 and WA92 data. This conclusion differs from the one of ref. [60] for the WA92 result. This 
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Figure 15: Azimuthal correlation for charm production in ttN collisions: NLO calcula- 
tion versus the WA75 and WA92 data. 



is because in ref. |58] the WA92 collaboration has improved the study of correlations with respect 



to ref. |n[ by considering a wider set of correlation variables and by improving the statistics 
by a factor of 5. WA92 and WA75 data now appear to be consistent with each other. As is 
apparent from fig. [15], the acceptable value of (A; 2 ) = 1 GeV 2 is required to describe the data. 
This is encouragingly consistent with the results we obtained form the study of the inclusive-pr 
distributions! 

The WA75 collaboration, and recently the WA92 collaboration, published in refs. [[K3|, Q the 
distribution of the transverse momentum of the heavy-quark pair. This is yet another indepen- 
dent variable, which is sensitive to both intrinsic k T and fragmentation effects. The theoretical 
prediction supplemented with a fc T -kick with (A; 2 ,) = 1 GeV 2 cannot reproduce the WA75 data. 
This is a typical example of why, until recently, no coherent overall picture of the mechanisms at 
play in charm production was available. The most recent results from WA92, however, confirm the 
indications of the QCD calculations supplemented with fc T -kick and fragmentation, as displayed 
in fig. |TB[ Similar conclusions can be drawn from the study of photoproduction correlation data 
(the A0 distributions are shown in fig. pj], the distribution in the transverse momentum of the 
heavy-quark pair in fig. |T^). In all cases (A; 2 ) values of the order of 1 GeV 2 are necessary (and 
sufficient!) to describe the data. 

Conclusions: 

• For the first time since the beginning of charm production studies we have reached a com- 
pelling and coherent picture of the production dynamics in fixed-target experiments. Large 
data-sets and state-of-the-art detector technology have been fundamental to make outstand- 
ing progress in the quality of the data! 

• The production dynamics are well described by NLO QCD, supplemented with simple phe- 
nomenological descriptions of the NP effects: 

— fragmentation at large Pt 
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Figure 16: iVXO QCD result for the p%(QQ) supplemented with an intrinsic transverse 
momentum for the incoming partons, compared with the WA75 (left) and WA92 (right) 
data. ,. 
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Figure 17: Azimuthal correlation of DD pair versus the perturbative result in photopro- 
duction for the E687 (left) and NA14/2 (right) experiments. 

- intrinsic (h?) = 1 GeV 2 

• The parametrization of these effects satisfies the expected universality properties: no depen- 
dence on 

— beam type and energy 

- observable 

• Additional fine tuning will be possible when more data will be available. We can then expect 
to sort out the interplay between QCD and NP parameters in a more accurate fashion. E.g.: 

— larger a s — * smaller (kj) (from A<p spectrum) 
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Figure 18: Transverse momentum distribution of the DD pair versus the perturbative 
result for the E687 experiment. 

— smaller (kj) — > softer fragmentation (from spectra) 



3.2 Heavy-quark production at hadron colliders 

Hadron collisions are by far the largest source of heavy quarks available today. While the en- 
vironment of high-energy hadronic interactions does not allow to trigger on the largest fraction 
of charm and bottom produced, the production rates are so large that the number of recorded 
events allows today 6-physics studies that are competitive with those of e + e~ experiments. The 
introduction of new experimental techniques, such as the use of silicon vertex detectors, which 
enable the tagging of events containing bottom quarks |§, led in the recent years to high-statistics 
and low-background measurements of the 6-production properties over a large range of transverse 
momenta. 

From the point of view of QCD studies, heavy-flavour production in high-energy hadronic 
collisions has better potentials than in fixed-target experiments. The b quarks produced at large 
p T can be studied in perturbative QCD with smaller contamination from non-perturbative effects. 
For example, the impact of the initial-parton transverse momentum is largely reduced with respect 
to fixed-target charm production. Furthermore, the fragmentation properties of heavy flavours in 
high transverse momentum jets can be directly studied, since the transverse momenta involved 
are typically perturbative. 

In spite of all of this, the comparison of theory with data collected by the high-energy collider 
experiments (UAl, CDF and DO) shows some puzzling discrepancies, which will be reviewed in 
the rest of this lecture. In particular, we will concentrate on the inclusive pt distributions in the 
central and forward rapidity regions, and on the azimuthal correlations between b and b. 
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3.2.1 Comparison of data and theory 



The distribution most commonly studied by the hadron-collider experiments is the 6-quark differ- 
ential pt spectrum, integrated within a fixed rapidity range and above a given pr threshold (p™ n 

): 

o(pT)= I dy f . dpT—^-. (38) 
J\y\<y max Jp T >pT n dydp T 

The UAl experiment at the CERN pp collider used y ma x = 1.5, while the CDF and DO experiments 
at the Tevatron use y max = 1. 

The results obtained by CDF [B2| and DO p3] are shown in figs. compared with the 



theoretical predictions. For an easier evaluation of the results, we present them on a linear scale 
in the form Data/Theory, and include the UAl data [64| as well. We divided the data by our 
central theoretical prediction. The dot-dashed straight lines are constant fits to the ratios, weighed 
by the inverse of the experimental uncertainties. The upper and lower solid lines correspond to the 
upper and lower theoretical predictions divided by the central prediction. The lower prediction 
is obtained fixing the renormalization and factorizations scales to be equal to 2/x (with /j,q = 
Pt 2 + m 2 .), and using the MRSA' parton density set p5 |. The value of A^ s for this set (152 MeV) 
yields a value of a s significantly lower than that extracted from different observables ||66|| . To 
estimate the upper theoretical curve we use the set MRS125 |3"I|], for which a s (M z ) = 0.125, 



and a renormalizatoin scale equal to /io/2. Given the large values of pt probed by the collider 
experiments, the 6-mass dependence of the theoretical result is small. We chose as a default the 
value rrib = 4.75 GeV, for all our curves. 

The first important thing to notice in fig. [19] is that, independently of the input parameters 
chosen, the shape of the theoretical curves agrees well with the data. Secondly, the results at \fS = 
630 GeV are by and large consistent with those at yS = 1800 GeV. The average ratio Data/Theory 
measured by UAl differs by less than 10% from the ratio measured by DO, independently of the 
input parameters chosen. The difference is however between 40 and 50% if one uses the CDF 
data. If we naively average the Tevatron data, we conclude that the relative discrepancy between 
theory and data when changing the value of VS from 630 to 1800 GeV is about 30%, a number 
of the same order of the estimate of log(l/x) effects discussed in the first lecture. 

Independently of the beam energy, the data are higher by a factor of about 2 than the default 
prediction based on fi = fi . They are, however, well reproduced by the upper theoretical curve. 
Therefore, while the overall uncertainty of the theoretical prediction due to the scale choice is 
large, there is currently no serious inconsistency between theory and data for the inclusive pt 
distribution of b quarks at the Tevatron. The 30% discrepancy between the results of CDF and 
DO is comparable to the discrepancy between the extrapolation of the UAl data to CDF, while 
UAl and DO data agree at the level of 10%. 

Nevertheless, the discrepancy between CDF and DO remains quite puzzling, as it is not con- 
sistent with the quoted systematic and statistical errors. In addition to this puzzle, a recent 
measurement of forward b production by DO |67j (performed using muons in the rapidity range 
2.4 < |?7 M | < 3.2) indicates a rate in excess by a factor of approximately 4 over the central theoret- 
ical prediction. The results are shown in fig. ^0[ Certainly no phenomenon of perturbative origin 
(such as higher-order corrections) could explain such a discrepancy. 

Similar inconsistencies are found in the study of correlations. A measurement has been recently 
reported by CDF, using muons plus tagged b jets |j6"8| . The shape of the A0 distribution is in 
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Figure 19: Linear comparison between experimental data and theory for the integrated 
b-quark pt distribution. See the text for details. 

reasonable agreement with QCD, while the distributions of the jet E T and of the muon are not 



see fig. 21). On the other hand, both CDF and DO recently presented studies based on samples 



of high- mass dimuons [|69"1 , fTOfl . The shapes for both and p\ for a given p b T are in this case well 



reproduced by theory, within the large uncertainties (fig. [22]) . 

Contrary to the case of cc correlations measured in fixed-target experiments, the measurement 
at hadron colliders is sensitive to the modelling of the heavy-quark fragmentation, because of 
possible trigger biases. A harder (softer) fragmentation function would enhance (decrease) the 
efficiency for the detection of the softest of the two b quarks. These effects could have an impact on 
the distributions reported by CDF and DO. CDF explored the effects of changes in the e& parameter 
within the Peterson fragmentation model, finding them negligible. As we argued earlier, it cannot 
be excluded that a systematic study of other possible parametrizations for the fragmentation 
modelling could lead to significant effects. Also the possible effects of the k T kick have been studied 
by CDF |p8| , with the conclusion that not even an average k T as large as 4 GeV could improve 
the agreement between theory and data in the case of the muon plus 6-tagged jet measurement. 

In summary, while the overall features of b production at the Tevatron are well described by 
NLO QCD, there are still some apparent inconsistencies which go beyond the quoted theoretical 
and experimental uncertainties. I will finish these lectures by presenting some studies which 
have been reported recently, and others which are in progress, aimed at a better assessment of the 
theoretical systematics. I will first present studies related to the resummation of large perturbative 
logarithms, and then some considerations on the effects of non-perturbative fragmentation. 
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Figure 20: Comparison between DO data and theory for the pp — > (b —* fi) + X cross- 
section at large rapidity. 

3.2.2 Resummation of large-p^ logarithms 

Two different groups have approached this problem in the recent past. Cacciari and Greco p0| 
have folded the NLO cross-section for production of a massless parton i (i = g, q) J7TJ with the 



NLO fragmentation function for the transition i — > b [I6[. The evolution of the fragmentation 
functions resums all terms of order a s n log™ (pr/m) and a s n+1 log™ {px/m). All the dependence on 
the 6-quark mass lies in the boundary conditions for the fragmentation functions. This approach 
ensures a full NLO accuracy, up to corrections of order m 2 / (m 2 +p 2 )- In particular, this formalism 
describes NLO corrections to the gluon splitting process, which in the O(af) calculation is only 
included at the leading order. One can verify, by looking at the Born-level production process as 
a function of the quark mass, that, in order for the mass corrections not to exceed the 10% level, 
it is however necessary to limit the applications of this formalism to the region of pt > 20 GeV. 

Figure [2B| shows the differential 6-quark p? distribution obtained in the fragmentation-function 
approach, compared to the standard fixed-order NLO result. Several features of this figure should 
be noticed. To start with, the scale dependence is significantly reduced with respect to the fixed- 
order calculation. Furthermore, in the range of applicability of this calculation (i.e. pt > 20 GeV) 
the result of the fragmentation-function approach lies on the upper side of the fixed-order NLO 
calculation. The resummed calculation is, however, always within the uncertainty band of the 
fixed-order one. Finally, notice that the overall effect of the inclusion of higher-order logarithms 
is a steepening of the pt spectrum, as is natural to expect. 



Another calculation has recently appeared, by Scalise, Olness and Tung [21]. In this approach 



the authors employ a strategy developed in the case of DIS in refs. [ 72 , [73]] . Here initial- and final- 
state mass singularities are resummed as in the fragmentation-function approach, and the result 
is then matched in the low-p-r region to the fixed-order NLO massive calculation. At large pt this 
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Figure 21: CDF results on the bb correlations using /i+b-tagged jet final states jpS\j, 
compared to NLO QCD. Solid points correspond to the default theoretical prediction, with 
scales fj, = empty circles correspond to the difference between the choice [i = Ho/2 and 
H = [j,q. Top figure: tagged-jet E T distribution. Central: muon px spectrum. Bottom: 
azimuthal correlations. 

calculation does not include as yet, however, the full set of NLO corrections to the hard-process 
matrix elements. 



The preliminary numerical results of this study |2T| are consistent with those of the approach 
by Cacciari and Greco; in particular, they support the conclusion that in the range explored by 
the Tevatron experiments the resummed results are consistent with the fixed-order ones, provided 
a scale fi of the order of yUo/2 or smaller is selected. 



3.2.3 Non-perturbative fragmentation 

We start with a remark: the use of the Peterson fragmentation function might not be justified in 
the context of hadroproduction of heavy quarks. As a simple observation, we point out here that 
while the measurement of heavy-quark spectra in e + e~ data is mostly sensitive to the first moment 
of the fragmentation function, corresponding to the average jet energy carried by the heavy hadron, 
the pt distributions in hadronic collisions are sensitive to higher moments of the non-perturbative 
fragmentation function. In fact, assuming for simplicity a perturbative pt spectrum of the form: 

^ = — , (39) 

dpT PT n 
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Figure 22: Results on bb correlations using dimuon final states, compared to NLO QCD. 
Left: azimuthal- correlation results from DO. Right: px distribution of the slowest muon 
from CDF. 

it is easy to prove that the resulting hadron spectrum, after convolution with a given fragmentation 
function f(z), will be given by: 

= — x / dzz n ~ l f{z). (40) 
dp T p T n J w K J 

In the case of the Tevatron, n turns out to be approximately 5. It is not unlikely, therefore, that 
alternative models for the non-perturbative fragmentation of heavy quarks, which could equally 
well fit the e + e~ data, could give rise to significant differences when applied to production in 
hadronic collisions. We also remark that the gluon component of the fragmentation function is 
not important in e + e~ physics, while it may be crucial in hadroproduction. 

Another important point concerns the correlation which exists between the chosen value of a s , 
the accuracy at which the perturbative fragmentation function is evaluated, and the value of the 
parameters describing the non-perturbative fragmentation extracted from the comparison with 
the e + e~ data. A larger value of a s , and a NLO calculation, give rise to a softer perturbative 
fragmentation function than would be obtained from a smaller value of a s or a LO accuracy. As 
a result, the non-perturbative fragmentation function will be harder. This was already noticed in 
the first lecture, where we compared the values of the Peterson parameter e extracted from a LO 
and a NLO fit to the Argus and Opal D-production data (e c = 0.06 at LO, e c = 0.01 at NLO). 

A complete study of the effect of these differences for b production at the Tevatron has not 
been carried out as yet. As a crude estimate of the effects, however, we can just compare the pr 
spectra obtained by convoluting the fixed-order NLO calculation of the b-px distribution with NP 
fragmentation functions fitted at the LO and at the NLO. In fig. |24| I plot the ratio of the 5-meson 
Pt distribution obtained by folding with the Colangelo-Nason fragmentation function (with the 
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Figure 23: Differential b-quark distribution: comparison of the fixed-order NLO 
calculation (NDE stands for with the fragmentation-function approach tfjQj ; fio is the 
scale at which the boundary conditions for the fragmentation functions are set. 

parameters a and (3 fitted at NLO), relative to that obtained by using the Peterson fragmentation 
function and e c = 0.06, the default choice of CDF and DO. This is done for two different rapidity 
intervals, one describing central production, the other describing forward production. 

In the central region we obtain differences of up to 30%, which are consistent with the excess 
in rate reported by CDF. In the forward region the effect is larger, up to a factor of 50%, but not 
large enough to explain the deviation seen by DO using their forward-muon data. 

We summarize the main conclusions of the studies presented in this section: 

1. There is good agreement between the shape of the 6-quark p T distribution predicted by NLO 
QCD and that observed in the data for central rapidities. 

2. Although the data are higher by a factor of approximately 2 with respect to the theoretical 
prediction with the default choice of parameters, extreme (although acceptable) choices of 
Af s and of renormalization and factorization scales bring the theory in perfect agreement 
with the data of UA1 and DO, and within 30% of the CDF measurements. 

3. The choice of low values of /zr and /ip, which helps the agreement between theory and data, 
is favoured by studies of higher-order logarithmic corrections. 
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4. The 30% difference between CDF and DO cannot be explained by theory, and will need to 
be understood before further progress can be made. 

5. Forward production of b quarks at DO indicates a larger discrepancy between theory and 
data, which cannot be accomodated by higher-order perturbative corrections, nor by the 
systematic uncertainties induced by the choice of non-perturbative fragmentation functions. 

As should be clear from the list above, more work is necessary both on the experimental and on 
the theoretical sides to improve our understanding of bottom production at high energy. To which 
extent the current inconsistencies will affect future physics plans at the Tevatron and at the LHC 
it is not clear. It is unlikely however that measurements such as CP asymmetries will be affected 
by these issues, since they rely on ratios of rates. Nevertheless the study of b production remains 
an important test of our ability to perform solid perturbative predictions, and any progress which 
can be made in the coming future will certainly benefit our understanding of QCD at high energies. 
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